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ABSTRACT 


Detection of primary brain tumors is inspired by the necessity of high accuracy 
as it deals with human life. Various imaging modalities techniques have 
incarnated as a tool in diagnosis and treatment domain. Yet, experienced and 
competent medical practitioners for the proper interpretation are still required. 
Thus, the involvement of information technology is highly demanded in 
introducing reliable and accurate computer systems. This study presents an 
algorithm for nosologic segmentation of primary brain tumors on Magnetic 
Resonance Imaging (MRI) brain images. Nosologic refers to the classification 
of diseases that can facilitate the diagnosis of neurological diseases. The 
purpose of segmentation is to highlight the tumor areas, whereas classification 
is used to identify the type of the primary brain tumors. For this purpose, an 
algorithm which hybridized the Grey Level Co- occurrence Matrices (GLCM), 
Intensity Based Analysis (IBA), Adaptive Network-based Fuzzy Inference 
System (ANFIS) and Particle Swarm Optimization (PSO) Clustering 
Algorithm (CAPSOCA) is proposed. The combination of several computer 
vision techniques is aim to deliver reproducible nosologic segmentation of 
primary brain tumors which are gliomas and meningiomas. The performance 
of the CAPSOCA is quantified by two measurements which are segmentation 
and classification accuracy. The segmentation accuracy is evaluated using 
comparison with ground truth approach. On the other hand, the classification 
accuracy is quantified using a truth table by comparing the classification 
outcomes with histopathology diagnosis. Upon the testing conducted, the 
CAPSOCA was proven to be an effective algorithm for nosologic 
segmentation of primary brain tumors. It appeared to return 88.09% of overall 
mean accuracy for gliomas segmentation, 86.92% of overall mean accuracy 
for meningiomas segmentation. In another note, 83.72% and 85.19% of 
classification accuracy for gliomas and meningiomas were observed. 
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1, INTRODUCTION 

A rapid development of medical imaging modalities such as X-ray, Computed Tomography (CT) 
scan, Ultrasound and Magnetic Resonance Imaging (MRI) has enabled the investigations of domains that 
encompass areas that were out of reaches for the human’s eyes. These modalities make it possible to explain 
the structures of organs and cells, enable the observation of the way they function, allows abnormalities or 
dysfunction detection as well as assisting in pathology diagnosis [1]. 
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The brain is defined as one of the most complicated, least accessible and prone to complex 
abnormalities in human organs [2]. Accordingly, medical imaging modalities are anticipated to contribute the 
primary benefits in overcoming these circumstances. The brain complexity can be expressed at a variety of 
scales. Therefore, a deeper understanding of the brain anatomical structures plays a crucial role in the search 
for more efficient brain lesions and diseases (brain tumor, infarction, haemorrhage, stroke and many other brain 
abnormalities). Detection of tumors in brain tissue area from different medical images is inspired by the 
necessity of high accuracy as it deals with human life [3]. Presently, various imaging modalities techniques 
have incarnated as a tool for the doctors and radiologists to help them in diagnosis and treatment domain. While 
these are highly accurate and fast, they still require experienced and competent medical practitioners for the 
proper interpretation [4]. Thus, computer aided diagnosis is significantly useful due to the fact that it could 
enhance the results of humans in such domain [5, 6]. 

The involvements of information technology in the development of various applications have 
completely changed the world. The noticeable reasons for the introduction of computer systems are reliability, 
accuracy, simplicity and ease of use [5-7]. Moreover, it has become almost compulsory to use computers to 
assist radiological experts in clinical diagnosis and treatment planning due to the increasing use of imaging for 
diagnosis, treatment planning and clinical studies [8, 9]. In the past several years, computer aided system has 
widely applied for many ranges of application such as a brain tumour [10], lung cancer [11] and breast cancer 
[12]. Computer vision is a branch of Artificial Intelligence (AI) that provides computers with the functions 
characteristic of human vision. Currently, computer vision has produced important applications in image 
processing perspective such as segmentation [13-15], feature extraction [16], texture classification [17] and 
many other applications. 

PSO has been observed to be successfully applied for clustering approach in many research and 
application areas, as well as various representations of methods for image segmentation. It has also gained 
increasing popularity in recent years due to its ability to solve efficiently and effectively an increasing number 
of problems in science and engineering as in these numerous studies [18-20]. It is reported that the PSO 
outperformed other methods such as Fuzzy c-Means (FCM) [21], K-means, Genetic Algorithm (GA) and 
Differential Evolution (DE) [22], and even the other Swarm Intelligence (SI) algorithms which are Ant Colony 
Optimization (ACO), Artificial Bee Colony Algorithm (ABCO) and Bacterial Foraging Algorithm (BFOA) 
[20]. Therefore, this paper presents a hybrid algorithm known as CAPSOCA for nosologic segmentation of 
primary brain tumors on Magnetic Resonance Imaging (MRI) brain images. The algorithm is designed as a 
PSO-based clustering algorithm which combines several computer vision techniques which are GLCM, IBA, 
and ANFIS. The GLCM and IBA acted as feature extraction techniques in determining the patterns and 
characteristics for the various types of brain tissues. The basic concept is that the local textures in the images 
can reveal the biological structures of gliomas and meningiomas. Alternatively, the ANFIS is implemented to 
classify the type of primary brain tumors. The encapsulation of all the implemented techniques into the PSO- 
based clustering algorithm produced a hybrid technique called as CAPSOCA. The organization of the rest of 
this paper is as follows: Section 2 presents our proposed research method, including the brief description on 
data collection, and the CAPSOCA structure. Section 3 presents our results and discussions. Finally, we present 
our conclusions in Section 4. 


2. RESEARCH METHOD 

The data collection involved 167 Fluid Attenuated Inversion Recovery (FLAIR)-MRI brain images 
of gliomas and menigiomas which were collected from Hospital Sungai Buloh, Selangor. The aim of this 
study is to generate an innovative non-invasive method for nosologic segmentation of primary brain tumors 
in the MRI field. The concept is mainly focused on the development of a comprehensive algorithm which 
combines several computer vision techniques for nosologic segmentation of primary brain tumors which are 
gliomas and meningiomas. The basic concepts of the PSO, the proposed hybrid technique of CAPSOCA, 
including the description of all techniques involved are explained in the next subsection. 


2.1. The PSO 

PSO is an evolutionary computation technique proposed by Kennedy and Eberhart in 1995. The basic 
idea of PSO 1s inspired by social behaviour of bird flocking, fish schooling and swarm theory. It is a SI based 
algorithms that have been applied in numerical optimization problems [23]. PSO is a very popular SI algorithm 
for global optimization over continuous search spaces. It has attracted the attention of several researchers which 
resulting in a huge number of variants of the basic algorithm as well as many parameter automation strategies 
ever since its advent in 1995 In the past several years, PSO has been successfully applied in many research and 
application areas. One of the advantages of PSO is that it 1s easier to implement and there are very few 
parameters to adjust [24]. It has also demonstrated as a faster and cheaper way of obtaining better results as 
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compared with other methods [25]. Reference [26] proposed a particle swarm based segmentation algorithm 
for automatically grouping the pixels of an image into different homogeneous regions. They used PSO for the 
fuzzy segmentation of images since the pixels in an image generally have high correlation and have proven 
that PSO performed well. 

In implementing the PSO, several matters should first be considered in order to apply the PSO idea 
such as the representation of initial population, representation of position and velocity strategies, fitness 
function identification and the limitation. PSO is defined by considering five essential parameters as tabulated 
in Table 1. 


Table 1. PSO parameters 


Parameters Description 
Particle candidate solution to a problem 
Velocity the rate of position change 
Fitness the best solution achieved 
pceurrent a position value of the current particle 
pbest present best, the best value obtained in previous particle 
gbest global best, the best value obtained so far by any particle in the population 


2.2. The CAPSOCA 

The CAPSOCA is proposed in which several computer vision techniques are combined throughout 
the PSO-based algorithm. The combination of GLCM and IBA is used to analyze the texture of the brain 
component categories. Whereas, the ANFIS plays the texture classification roles which classify the texture to 
the category it belongs to. Finally is the implementation of PSO algorithm which acted as the search and 
optimization method in finding the best pixels for segmentation, and also as the outline for the whole nosologic 
segmentation algorithm. The pictorial representation of the proposed nosologic CAPSOCA segmentation 
model is illustrated in Figure 1. 








Figure 1. Proposed nosologic capsoca segmentation model 


The six primary steps of CAPSOCA segmentation are elaborated as below: 
a. Step 1: Generation of Particles 

PSO 1s initialized with a swarm of random solutions called particles. Each particle maintains its own 
position, velocity and its current best position. The algorithm is initialized to contain 500 points of particles 
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with random position and velocity. The points had been randomly selected in the X-axis value within the image 
width while the Y-axis value within the image height. Therefore, the computer will automatically generate 500 
points from the image randomly without considering the position of the points. 
b. Step 2: GLCM and IBA Fitness Examination 

After the swarm of particles is generated, each particle’s fitness needs to be examined based on the 
Region of Interest (ROI) table which acted as the objective function. An objective function refers to a particular 
type of function that quantifies the solution of an optimization problem. This function is specific to each 
optimization problem and its suitability is essential to contribute a success and performance of the particular 
optimization process. An ideal objective function correlates closely with the goal of a particular algorithm. In 
some cases, it is very hard or impossible to come up even with a guess of what fitness function definition might 
be. In this study, the ROI table combines both GLCM and IBA parameters. The pattern and characteristics of 
each particle and its 3x3 matrix surrounding pixels are examined according to the parameters values defined 
in the reference table. The objective function used is formulated as in (1): 


Min (Obj Funcpc) 
Max (ObjFuncpc) 
Mean (ObjFuncpc) (1) 


where p = 0, 1, 2, 3, 4, 5, 6 and c = 0, 1, 2, 3. p represents the seven parameters of GLCM and IBA which are 
energy, entropy, contrast, dissimilarity, min, max and mean. Whereas, the c represents the four brain 
component categories of dark gliomas, light gliomas, dark meningiomas and light meningiomas. 
The parameters are further described as in Table 2. 


Table 2. Parameters of reference table 


Method Parameter Description 
Energy Provides the sum of squared elements in the GLCM 
Entropy Measures the uncertainty or disorder 
Contrast Measures the local variations in the GLCM 
GLCM Dissimilarity Measures the discrepancy between the two objects based on several 


features. Dissimilarity may also be viewed as a measure of disorder 
between two objects 


Min Minimum grey level pixel value occurred 
IBA Max Maximum grey level pixel value occurred 
Mean Mean of grey level pixel value occurred 


These values are then matched with the same parameters of the current particle so that the current 
particle could be classified into the particular brain component it belongs to. The classification of each particle 
is performed using ANFIS technique. 

c. Step 3: ANFIS Classification 

The classification of each fit particle is achieved using a classification method of ANFIS. The step by 
step for each layer of ANFIS classification is presented in a flow chart as shown in Figure 2. After a particular 
fit particle 1s classified by the ANFIS, the particle is then segmented according to the brain component 
categories it belongs to. Each region of tumors is assigned to different colors with the purpose that the 
segmentation of each region is noticeable. 

d. Tumor Segmentation 

Segmentation of medical images holds an important position in the area of image processing. It is 
used to extract information from complex medical images and it has wide application in the medical field. The 
main objective of tumors segmentation in this study is to segment the region of tumors which are gliomas and 
meningiomas. Table 3 tabulates a sample of CAPSOCA segmentation for gliomas and meningiomas. 

e. Step 5: Position and Velocity Update 

Next, each current pbest particle is evaluated in the same way as the particles generated initially. Each 
best fitted particle among all the particles is stored as global best (gbest). This process is done iteratively which 
the fitness’s, position and velocity are updated in every iteration if necessary. The iteration of update rules of 
position and velocity leads to the exploration of the whole regions that turn out to be the final outcomes. 
Velocity refers to the rate of position change which the particle’s position is depending on the velocity. In the 
first iteration, four neighbouring pixels of each particular pbest particle are considered. There are four 
coordinates of neighbouring pixels considered are (x, y-1), (x-1, y), (x, y+1) and (x+l, y). 
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Figure 2. Flowchart of ANFIS classification 


Table 3. Samples of CAPSOCA segmentation 
Category Original MRI Image CAPSOCA Segmentation 


Gliomas 


Meningiomas 
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f. Step 6: Determine Termination Criteria: 

The processes of fitness examination, velocity and position update are performed constantly until the 
termination criteria are met. These criteria are determined based on the problem to be optimized. Therefore, 
the process will halt and return the result when there are no unprocessed pbest and gbest pixels in the maximum 
velocity regions. Otherwise, the process will continue to start the fitness examination, velocity and position 
updating processes for the next ith particles. The maximum velocity is used to keep the particles from moving 
too far beyond the search space. The maximum velocity in this study is limited to the whole region of each 
MRI brain image. It is chosen for facilitating global exploration of particle’s position since too low maximum 
velocity region might lead to the difficulties of particles in exploring the optimal regions. 


2.3. Evaluation of Accuracy 

The CAPSOCA segmentation quality and the level of accuracy for gliomas and meningiomas 
are quantified using statistical values which are the mean percentage of sensitivity, specificity and accuracy as 
in (2)-(4). 


Sensitivity = ee ee, 100% (2) 


True Positive+False Negative 


True Negative 


Specificity = X 100% (3) 


True Negative+False Positive 


True Positive+True Negative 


Accuracy = X 100% (4) 


True Positive+True Negative+False Positive+False Negative 


On the other hand, the primary brain tumors classification accuracy 1s evaluated using a truth table. A 
truth table is a useful assessing method for determining the validity or invalidity of any form of argument. 
Pixel-based diagnosis 1s proposed to classify the types of primary brain tumors based on the tumor- segmented 
pixels. It 1s done by comparing the percentage values of pixel-based diagnosis for both gliomas and 
meningiomas. The percentage value of pixel-based is calculated by the numbers of tumor-segmented pixels 
which overlapped with the ground truth. The pixel-based diagnosis is determined according to the highest 
percentage of overlapped pixel value between both types of primary tumors. Table 4 tabulates a number of 
samples on how the pixel-based diagnosis 1s identified. 


Table 4. Samples of pixel-based diagnosis 


CAPSOCA Ground Truth % of Overlapped Pixels *Pixel-based Histopathology Level of 
Segmentation Gliomas Meningiomas Diagnosis Diagnosis Confidence 
43.6 89.4 Meningiomas Meningiomas TRUE 
84.4 = pe Gliomas Gliomas TRUE 





*Pixel-based diagnosis is identified based on the highest value of the percentage of overlapped pixels 


Next, the level of confidence measured the accuracy of the pixel-based diagnosis by comparing the 
outcomes with histopathology diagnosis collected from the patient’s record. The level of confidence 1s divided 
into two conditions which are TRUE or FALSE. Consequently, the percentage of accuracy for the pixel-based 
diagnosis model is calculated according to the number of TRUE pixel-based diagnosis to the total number of 
images as in (5). 


No.of TRUE Pixel—based Diagnosis 


% of Classification Accuracy = Naaru 


(5) 
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3. RESULTS AND ANALYSIS 

The performance of the CAPSOCA is quantified by two measurements which are segmentation 
accuracy and classification accuracy. The segmentation accuracy is evaluated using comparison with ground 
truth approach. Table 5 tabulates the CAPSOCA segmentation summary for both gliomas and meningiomas. 


Table 5. Summary of gliomas and meningiomas segmentation 


Category Basic Quantifier 
Sensitivity % Specificity % Accuracy % 
Gliomas 78.55 88.62 88.09 
Meningiomas Pee 88.34 85.92 


As seen from the Table 5, the CAPSOCA is observed to produce good segmentation results for 
gliomas. The statistic shows that the sensitivity, specificity and accuracy mean values produced consistent 
outcomes which are 78.55%, 88.62%, and 88.09% respectively. These proved that the CAPSOCA 1s capable 
of producing good quality of segmentation in gliomas. 

On the other hand, the CAPSOCA segmentation outcomes for meningiomas are also cannot be 
underestimated since it produced high mean values of specificity and accuracy, which are 88.34% and 85.92%. 
However, a slightly low mean value of sensitivity is observed. This makes the segmentation of meningiomas 
is not as good as gliomas. The mean of sensitivity value proved that a high rate of False Negative might occur 
in segmentation which the CAPSOCA is incapable to segment the area it supposed to be. This is caused by 
confusion in distinguishing the texture homogeneity between the meningiomas with other anatomical brain 
structure such as the membrane area. 

On the other hand, the classification accuracy is quantified using a truth table by comparing the 
classification outcomes with histopathology diagnosis collected from the patient’s record. Table 6 
summarizes the truth table results for gliomas and meningiomas classification. 


Table 6. Summary of gliomas and meningiomas classification 
Category No.of No.of TRUE No. of FALSE % of 
Images Classification Classification | Accuracy 
Gliomas 86 72 14 83.72 
Meningiomas 81 69 12 85.19 


Referring to Table 6, the primary brain tumors classification accuracy is evaluated to 86 images of 
gliomas and 81 images of meningiomas. The percentage of accuracy for meningiomas is observed to produce 
higher percentage value of 85.19% as compared to the gliomas which returned a slightly lower value of 83.72% 
of accuracy. Based on the performance, it can be concluded that the CAPSOCA showed a great potential for 
primary brain tumors classification. 


4. CONCLUSION 

This paper presents a CAPSOCA algorithm for nosologic segmentation of primary brain tumors on 
Magnetic Resonance Imaging (MRI) brain images. The application to 167 Fluid Attenuated Inversion 
Recovery (FLAIR)-MRI brain images has been successful. The performance of the CAPSOCA towards the 
gliomas and meningiomas had been evaluated by two measurements which are segmentation accuracy and 
classification accuracy. The CAPSOCA is observed to produce good segmentation results for gliomas as 
compared to the meningiomas. The performance obtained exhibit a little variation in segmenting the 
meningiomas due to False Negative value which affect the mean of sensitivity value. In contrast, the 
CAPSOCA returned high classification accuracy for both gliomas and meningiomas. Therefore, it can be 
concluded that the proposed hybrid algorithm of CAPSOCA 1s found to be successful as it showed a great 
potential for nosologic segmentation of primary brain tumors. It is expected to offer anew way to produce high 
resolution nosologic segmentation of MRI brain images which may help medical practitioners in making non- 
invasive decision steps in the future. 
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